Low scale saturation of Effective NN Interactions and their Symmetries 
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The Skyrme force parameters can be uniquely determined by coarse graining the NN interactions 
at a characteristic momentum scale. We show how exact V[ ow k potentials to second order in momenta 
are accurately and universally saturated with physical NN scattering threshold parameters at CM 
momentum scales of about A = 250MeV for the S-waves and A = lOOMeV for the P-waves. The 
pattern of Wigner and Serber symmetries unveiled previously is also saturated at these scales. 
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The derivation of effective interactions from NN dy- 
namics has been a major task in Nuclear Physics 
ever since the pioneering works of Moshinsky [l| and 
Skyrme @. The use of those effective potentials, referred 
to as Skyrme forces, in mean field calculations can hardly 
be exaggerated due to the enormous simplifications that 
are implied as compared to the original many-body prob- 
lem |3H6|. Similar ideas advanced by Moszkowski and 
Scott [3_have become rather useful in Shell model calcu- 
lations @, II]. The Skyrme (pseudo)potential is usually 
written in coordinate space and contains delta functions 
and its derivatives Q. In momentum space it corre- 
sponds to a power expansion in the CM momenta (p' 
and p) corresponding to the initial and final state re- 
spectively. To second order in momenta the potential 
reads 

V(p',p) = J d 3 a;e- ix -( p '- p V(x) 
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where P a = (1 + <j\ ■ <y%)j1 is the spin exchange operator 
with P a = — 1 for spin singlet S — and P„ = 1 for spin 
triplet S = 1 states. In practice, these effective forces are 
parameterized in terms of a few constants which encode 
the relevant physical information and should be deduced 
directly from the elementary and underlying NN interac- 
tions. Unfortunately, there is a huge variety of Skyrme 
forces depending on the fitting strategy employed (see 
e.g. US, 03). This lack of uniqueness may indicate that 
the systematic and/or statistical uncertainties within the 
various schemes are not accounted for completely. In- 
terestingly, the natural units for those parameters have 
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been outlined in Ref. [12], |l3[ yiel ding the correct order of 
magnitude. A microscopic basis \l4 Il5j for the Density 
Functional Theory (DFT) approach has also been set up, 
but still uncertainties remain. 

Although the pseudo-potential in Eq. (TTJ) may be taken 
literally in mean field calculations, due to the finite exten- 
sion of the nucleus, its interpretation in the simplest two- 
body problem requires some regularization to give a pre- 
cise meaning to the Dirac delta interactions. The stan- 
dard view of a pseudo-potential (in the sense of Fermi) 
is that it corresponds to the potential which in the Born 
approximation yields the real part of the full scattering 
amplitude. This is a prescription which implements uni- 
tarity, but necessarily fails when the scattering length is 
unnaturally large as it is the case for NN interactions. 
On the contrary, the Wilsonian viewpoint corresponds 
to a coarse graining of the NN interaction to a certain 
energy scale. There are several schemes to coarse grain 
interactions in Nuclear Physics. The traditional way has 
been by using the oscillator shell model, where matrix 
elements of NN interactions are evaluated with oscillator- 
constants of about b — 1.4 — 2fm 9]. A modern way of 
coarse graining nuclear interactions is represented by the 
Viowk method [l6[ (for a review see [l7|) where all momen- 
tum scales above 2fm _1 are integrated out. The recent 
Euclidean Lattice Effective Field Theory (EFT) calcula- 
tions (for a review see e.g. [Hj]), although breaking rota- 
tional symmetry explicitly, provide a competitive scheme 
where coarse grained interactions allow ab initio calcula- 
tions combining the insight of EFT and Monte-Carlo lat- 
tice experience, with lattice spacings as large as a = 2fm. 
These length scales match the typical inter-particle dis- 
tance of nuclear matter d = ~ 2fm. Actually, the 
three approaches feature energy-, momentum- and con- 
figuration space coarse graining respectively and ignore 
explicit dynamical effects below distances ~ b ~ 1/A ~ a 
which advantageously sidestep the problems related to 
the hard core and confirm the modern view that ab ini- 
tio calculations are subjected to larger systematic uncer- 
tainties than assumed hitherto. Clearly, any computa- 
tional set up implementing the coarse graining philoso- 
phy yields by itself a unique definition of the effective 
interaction. However, there is no universal effective in- 
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teraction definition. For definiteness, we will follow here 
the Vjowk scheme to determine the effective parameters 
because within this framework some underlying old nu- 
clear symmetries, namely those implied by Wigner and 
Serber forces, are vividly displayed [19l422l |. 

In the present paper we want to show that in fact these 
parameters can uniquely be determined from known NN 
scattering threshold parameters by rather simple calcu- 
lations by just coarse grain the interaction over all wave- 
lengths larger than the typical ones occurring in finite 
nuclei. As we will show, this introduces a momentum 
scale A in the 9 effective parameters to, 1,2, £0,1,2 and 
tu,T,V which allow to connect the two body problem to 
the many body problem. Going beyond Eq. ([TJ requires 
further information than just two-body low energy scat- 
tering, in particular knowledge about three and four body 
forces and their scale dependence consistently inherited 
from their NN counterpart. The finite fcf situation rele- 
vant for heavy nuclei and nuclear matter involves mixing 
between operators with different particle number and, in 
principle, could be conveniently tackled with the method 
outlined in Ref. [23j where the lack of genuine medium 
effects is manifestly built in. 

For completeness, we review here the Viowk ap- 
proach [24j in a way that our points can be easily 
stated. The starting point is a given phenomenological 
NN potential, V, and usually denominated bare potential, 
whence the scattering amplitude or T matrix is obtained 
as the solution of the half-off shell Lippmann-Schwinger 
(LS) coupled channel equation in the CM system 



T^(fc',fc;fc 2 ) = V^(fc',fc) 
M N dqq 2 



E 
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V'(k',q)T l i d (q,k;k 2 )(2) 



where J is the total angular momentum and I, I' are or- 
bital angular momentum quantum numbers p,p',q are 
CM momenta and Mjy is the Nucleon mass. The unitary 
(coupled channel) S-matrix is obtained as usual 



S?, l (p) = 6 l ,, l -i^T l J l (p,p) 



(3) 



Using the matrix representation S J = (M J — il)(M J + 
with (M J )t = M J a hermitian coupled channel 
matrix, at low energies the effective range theory for cou- 
pled channels reads 

p l+l ' +1 Mi!,(p) = -(a-%, + ^(r^p 2 + (v^p* + . . .(4) 

which in the absence of mixing and using Si (p) — e 2lS ' ( p ' 
reduces to the well-known expression 



P 2l+1 cot 6i(p) = — + \r lP 2 + vip A + 
ai 2 



(5) 



An extensive study and determination of the low en- 
ergy parameters for all partial waves has been carried 



out in Ref. [25[ for both the Nijmll and the Reid93 po- 
tentials [2(1 | yielding similar numerical results. Dropping 
these coupled channel indices for simplicity the Vi ow k po- 
tential is then defined by the equation 

T(k',k;k 2 ) = V lowk (k',k) 



Mjv dq q 2 
(2tt) 3 k 2 - q 2 



Mowk(fc', q)T(q, fc; fc 2 ) , (6) 



where (fc, fc') < A. We use here a sharp three-dimensional 
cut-off A to separate between low and high momenta 
since our results are not sensitive to the specific form of 
the regularization. Thus, eliminating the T matrix we get 
the equation for the effective potential which evidently 
depends on the cut-off scale A and corresponds to the ef- 
fective interaction which nucleons see when all momenta 
higher then the momentum scale A are integrated out. It 
has been found [24j that high precision potential models, 
i.e. fitting the NN data to high accuracy incorporating 
One Pion Exchange (OPE) at large distances and de- 
scribing the deuteron form factors, collapse into a unique 
self-adjoint nonlocal potential for A ~ 400 — 450MeV. 
This is a not a unreasonable result since all the poten- 
tials provide a rather satisfactory description of elastic 
NN scattering data up to p ~ 400MeV. Note that this 
universality requires a marginal effect of off-shell ambi- 
guities (beyond OPE off-shellness) , which is a great ad- 
vantage as this is a traditional source for uncertainties in 
nuclear structure. Actually, in the extreme limit when 
A — > one is left with zero energy on shell scattering 
yielding T(fc, fc) - (2tt) 3 a Q / 'M N . 

Moreover, for sufficiently small A, the potential which 
comes out from eliminating high energy modes can be 
accurately represented as the sum of the truncated orig- 
inal potential and a polynomial in the momentum [27| . 
However, as discussed in [20l | a more convenient repre- 
sentation is to separate off all polynomial dependence 
explicitly from the original potential 



Mowk(fc', fc) = Vk N (fc', fc) + Vc A T (fc', fc) , 



(7) 



with (fc, fc') < A, so that if VQ T (k',k) contains up to 
0(p n ) then VWO', k) starts off at 0(p n+1 ) : i.e. the next 
higher order. This way the departures from a pure poly- 
nomial may be viewed as true and explicit effects due to 
the potential and more precisely from the logarithmic left 
cut located at CM momentum p = im/2 at the partial 
wave amplitude level due to particle exchange with mass 
m. Specifically, 



V(^p(fc',fc) = fc'fc''' C L f (A) + Dj' (A) (fc 2 



fc' : 



,(8) 



where the coefficients C l j (A) and D l j (A) include all con- 
tributions to the effective interaction at low energies. Al- 
though we cannot calculate them ab initio we may relate 
them to low energy scattering data, in harmony with the 
expectation that off-shell effects are marginal. Not sur- 
prisingly the physics encoding the effective interaction 
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in Eq. (JSJ) will be related to the threshold parameters 
defined by Eq. (@}. Thus, the relevance of specific mi- 
croscopic nuclear forces to the effective (coarse grained) 
forces has to do with the extent to which these threshold 
parameters are described by the underlying forces and 
not so much with their detailed structure. We will dis- 
cuss below the limitations to this universal pattern. 

Using the partial wave projection [28[ we get the po- 
tentials in different angular momentum channels. These 
parameters can be related to the spectroscopic notation 
used in Ref. [29[ . The S- and P-wave potentials are 
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The 9 effective parameters depend on the scale A and can 
be related to the effective force representation io,i,2 ; ^0,1,2 
and ty,u,S of Eq. (QJ by the following explicit relations, 
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The corresponding T-matrices are conveniently solved 
by factoring out the centrifugal terms 

r,/,,(fc', k) = k l k' v [tl^p) + u V j(p)(k 2 + k' 2 ) + . . . ](ii) 

which reduce the LS equation to a finite set of alge- 
braic equations which are analytically solvable (see e.g. 
Ref. [3l| and references therein). In the simplest case 
where only the C's are taken into account the explicit 
solutions for S- and P- waves are, 
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FIG. 1: (Color on-line) Counterterms for the S- (in MeVfm 3 , 
upper panel) and P-waves (MeVfm , lower panel) as a func- 
tion of the momentum scale A (in MeV). C's from Eqs. ([9]) 
solving Eqs. (((j} including the D's using just the low energy 
threshold parameters from Ref. |25l ] (Thick Solid) . C's ex- 
tracted from the diagonal Vi ow k(p, p) potentials [2J| at fixed 
A = 420MeV for the Argonne-V18 Q (dashed). C's for P- 
waves including D-terms without mixings (thick dotted). 



where ao(ai) is the scattering length (volume) defined 
by Eq. ©. The Eq. §FZ§ illustrates the difference be- 
tween a Fermi pseudo-potential and a coarse grained po- 
tential as the former corresponds to aoA <C 1 where 
Cs(A) ~ 16ir 2 ao/Miy. In the case aoA >• 1 one has 
instead Cs(A) ~ — 87r/(AfA). Full solutions including 
the D's are also analytical although a bit messier, so we 
do not display them explicitly. They rely on Eq. (Q} with 

ai So' a3 &_2« 3 Po> m Pn a3 P 2 ' ai Fu a E x , T3 Sl and n So 
(see Ref. [25| for numerical values for Nijmll and Reid93 
potentials). At the order considered here we just men- 
tion that while all P-waves constants run independently 
of each other with A the spin-singlet parameters Cig , 
Di 5 on the one hand and the spin-triplet parameters 
Cag , D3 Sl and Ce x on the other are intertwined. 

We now turn to our numerical results. As can be seen 
from Fig. [T] the comparison of contact interactions using 
threshold parameters with V[ ow k results evolved to A = 
420MeV [2J| (note the different normalization as ours) 
from the Argonne-V18 bare potential [30( are saturated 
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FIG. 2: (Color on-line) S — D waves mixing on C3 Sl (in 
MeVfm 3 ) as a function of the momentum space cut-off A (in 
MeV) in the 3 Si channel. We compare when including a) only 
C3 Sl (dashed), b) C3$ 1 and D3 Sl (dotted), c) C^^, D3 Sl and 
Dex (solid). Wigner symmetry is displayed when the mixing 
is included. See also Fig. [T] and Eq. (0 in the main text. 

for A = 250MeV for S-waves and for much lower cut-offs 
for P-waves. Note that this holds regardless on the details 
of the potential as we only need the low energy threshold 
parameters as determined e.g. in Ref. (25j |. The strong 
dependence observed at larger A values just reflects the 
inadequacy of the second order truncation in Eqs. (|9|). 
This also reflects in the 25% — 50% inaccuracy are off 
the exact Viowk of the D's themselves despite showing 
plateaus, and thus will not be discussed any further. 

The identity Ci So (A) = Cs Sl (A) for A > 250MeV fea- 
tures the appearance of Wigner symmetry as pointed 
out in Ref. but now we see that this does not de- 
pend on details of the force. Actually, the effect of the 
3 5i — 3 Diwave mixing represented by a non-vanishing off 
diagonal potential VeiOAp) becomes essential to achieve 
this identity (a fact disregarded in Ref. [IH). As can be 
seen from Fig. [2] there is a large mismatch at values of 
A ~ 200 — 300MeV when De ± is set to zero (and hence 
a Ex =0) as compared with the case De x 7^ 0. 

The scale dependence of the Skyrme interaction pa- 
rameters (not involving the D's) can be seen in Fig. [3] in 
comparison with the Viowk potentials [24| deduced from 
the Argonne-V18 bare potentials J30j. The plateaus ob- 
served in the different partial waves are corroborated here 
as well as a remarkable accuracy in reproducing the exact 
Vi ow k numbers. Moreover, the weak cut-off dependence 
of the spin orbit interaction observed in Fig. [3] suggests 
taking A — > in which case 

W ° = 2^ ( 2a3p ° + 3 " 3Pl ~ 5a3p ^ ' ( 13 ) 

which upon using Ref. [H[ yields W a = 72MeVfm 4 . This 
numerical value reproduces within less than 10% the ex- 
act Vi ow k value. As can be seen from Fig. [3] the effec- 
tive range correction r\ provides via additional D coeffi- 
cients the missing contribution. This is a bit lower than 



what it is found in phenomenological approaches from 
the p 3 /2 — Pi/2 level splitting in 16 [fj. In any case, the 
comparison with phenomenological approaches based on 
mean field calculations may be tricky since as already 
mentioned not all the terms are always kept, and selec- 
tive fits to finite nuclear properties may overemphasize 
the role played by specific terms. 

It has recently been argued that counterterms are fin- 
gerprints of long distance symmetries [l9l - l2l1 |. This re- 
markable result holds regardless on the nature of the 
forces and applies in particular to both Wigner and Ser- 
ber symmetries. We confirm that to great accuracy, Xq — 
(Wigner symmetry) and X2 = — 1 (Scrber symmetry). 
The astonishing large- N c (N c is the number of colours in 
QCD) relations discussed in Refs. [l9|-[22j provide a di- 
rect link to the underlying quark and gluon dynamics and 
after [34| suggests a 1/N% accuracy of the Wigner sym- 
metry in even-L partial waves. Wigner symmetry has 
proven crucial in Nuclear coarse lattice (a ~ 2fm) calcu- 
lations [18( in sidestepping the sign problem for fermions. 
As we see for the scales typically involved there this works 
with great accuracy already at A ~ 250MeV. Taking into 
account that we deal with low ener gies , it is thus puzzling 
that Chiral interactions to N 3 LO [33j having chiral cut- 
offs A x ~ 600MeV tend to violate Wigner symmetry in 
the Vi ow k sense, i.e. C* s ^ C* 5 , whereas smaller values 
A x ~ 450MeV are preferred. 

Within the low energy expansion we have neglected 
terms 0(p' 4 , p 4 , p' 2 p 2 ) which correspond to P-waves and 
S-wave range corrections. In configuration space this 
corresponds to a dimensional expansion, since 8{f 12) = 
0(A 3 ) and {P 2 ,£(f 12 )} - 0(A 5 ), {P\S(r 12 )} = 0(A 7 ). 
Within such a scheme going to higher orders requires 
also to include three-body interactions, ~ ^('"12)^(^13) = 
C(A 6 ). Actually, at the two body level there are more po- 
tential parameters than low energy threshold parameters. 
For instance, in the 1 Sq channel one has two independent 
hermitean operators, p' 4 + p 4 and 2p' 2 p 2 (which are on- 
shell equivalent), but only one vis threshold parameter 
in the low energy expansion (see Eq. (Q}). As it was 
shown in Ref. [35| (see also Ref. (36[) these two features 
are interrelated since this two body off-shell ambiguity 
is cancelled when a three body observable, like e.g. the 
triton binding energy, is fixed . An intriguing aspect of 
the present investigation is the modification induced by 
potential tails due to e.g. pion exchange which cannot be 
represented by a polynomial since particle exchange gen- 
erates a cut in the complex energy plane. The important 
issue, however, is that the low scale saturation unveiled in 
the present paper works accurately just to second order 
as long as the low energy parameters determined from 
on-shell scattering are properly reproduced. 
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